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Abstract We analyze a generic model of mesoscopic machines driven by the 
nonadiabatic variation of external parameters. We derive a formula for the 
probability current; as a consequence we obtain a no-pumping theorem for 
cyclic processes satisfying detailed balance and demonstrate that the rectifi- 
cation of current requires broken spatial symmetry. 
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Molecular motors and other mesoscopic machines perform their duties in 
the face of the violent thermal agitation of the microscopic world [3] . Living 
organisms employ a host of such machines to pump ions across cell mem- 
branes, transport cargo, and cause muscle contraction, among many other 
tasks [TTll 1 3 j . Remarkably, artificial nanoscale machines have in recent years 
been synthesized and investigated in laboratories around the world [T2ll21| . 
Apart from their intrinsic allure these machines suggest novel technologies 
such as microsopic particle segregators [TKHIE] - In view of the diversity of 
natural and artificial mesoscopic machines, there is growing interest in de- 
veloping theoretical frameworks to provide insight, guide experiments, and 
stimulate conceptual breakthroughs. 

The nanoscale machines found in biological systems are typically au- 
tonomous, powered by the consumption of chemical fuel (such as ATP) or 
some other steady supply of energy (e.g. sunlight). By contrast, many ar- 
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tificial machines are non- autonomous: external parameters, such as electric 
fields, temperature, pressure, and chemical potentials, are varied with time to 
produce the desired behavior. For such machines, a natural goal of theory is 
to describe how the system responds to these time-dependent perturbations. 

In this article we analyze stochastic pumps, a generic model of non- 
autonomous mesoscopic machines operating in the presence of random noise. 
We specifically focus on systems that evolve diffusively in one dimension - 
Brownian motors provide a concrete example [31IT1IMT7] - and we are inter- 
ested in the current J(x,t) that arises due to the "pumping" of external 
parameters A(i). (Discrete-state stochastic pumps have long been studied 
in biological systems [23]. For recent analyses of discrete-state stochastic 
pumps, see Refs. [T^J^O , 2 , 16 , 6J .) Parrondo [T5] has formulated this prob- 
lem similarly, and has solved for the current induced by adiabatic (i.e. slow) 
pumping. Here we analyze the nonadiabatic case in which the parameters are 
manipulated at arbitrary rates. We derive a nontrivial expression for J(x, t) 
(Eq.O and we explore its consequences. In particular, we obtain a geometric 
expression for the integrated current in the limit of slow pumping (Eq. 1 1 6|1 ; 
we extend a no-pumping theorem that we had previously obtained for dis- 
crete systems |16) ; and we provide a simple mathematical argument for the 
widely known fact that the generation of nonzero integrated current requires 
broken spatial symmetry [T71I9], 

Let us consider a system that can be modeled as a one-dimensional diffu- 
sion process [8] on the circle. The state of the process at time t is denoted by 
x(t) and takes values in the range x E [0,L], with the end points identified. 
This class of diffusion processes encompasses variables that are intrinsically 
periodic, such as the dihedral angle of a chemical bond, as well as extended 
reaction coordinates evolving in a periodic potential [17] • The probability 
density P(x, t) to observe state x at time t evolves according to the Fokker- 
Planck equation [18] 



-§- x ^) + ^B A (x) 



P(x,t) 



(1) 



= C A {x)P(x,t), 



with periodic boundary conditions. The dynamics depend on the external 
parameters A through the drift coefficient A A (x) and diffusion coefficient 
Ba{x) which are assumed to be periodic in x. The system is "pumped" by 
varying the parameters in time from t = to r according to a specified 
protocol A(t), from A(0) — el to A(r) = b, and we are interested in the 
flow of probability that arises as a result of this pumping. (For brevity, the 
time-dependence of A will be left implicit throughout the following.) 

Since probability is conserved, it is natural to cast Eq. [1] as a continuity 
equation, 

d d 

—P(x,t) + —J(x,t)=0, (2) 
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where the instantaneous current 



J(x,t) = A A (x) - ^-B A (x) 
= Ja{x)P{x,£) 



P(x,t) 

(3) 



is the rate of flow of probability, in the positive direction, past a fixed obser- 
vation point x at time t. We will also investigate the integrated current 

$(x) = [ dtJ(x,t), (4) 
Jo 

which measures the net flow of probability past the point x during the time 
interval < t < t. 

To analyze such diffusion processes, it is convenient to define two auxiliary 
functions iPa(x) and ip A (x), which are not necessarily periodic in x: 

f x A\(v) 

tp A {x)=lnB A (x)- dy-^-=]nB A (x)+7p A (x). (5) 
Jo B A{y) 

We will refer to ip A (x) as the potential, in view of the role this function plays 
when the dynamics satisfy detailed balance (see below). Observe that A A (x) 
and B A (x) can be reconstructed from <p A (x) and ^ A {x), that is Eq. [5] can 
be inverted. In other words the diffusion process is characterized equally well 
by the auxiliary functions as by the drift and diffusion coefficients. 

For every fixed A, there exists a unique stationary distribution P A (x), 
satisfying C A P A = 0, with stationary current J A = J A {x)P A {x) Q. In the 
special case that ^a(O) = (p A (L), it follows from Eqs. [I] and [3] that P A (x) oc 
e -tPA(x) an( j js^ _ o j n thjg cas6j we w i\\ sa y t na t £ A satisfies detailed balance 
and we will view the stationary distribution as an equilibrium distribution, 

For the setup we have just outlined, there are two distinct mechanisms 
to produce currents [5]. First, when detailed balance is broken the stationary 
distributions themselves support non-zero currents, J A ^ 0. Second, if the 
system is driven away from the stationary distribution by varying A with 
time, currents arise due to the resulting redistribution of probability. In this 
article we derive an explicit decomposition of the total current into these two 
contributions: 

J(x,t) = J s A + [ dx'V A (x,x')P(x',t), (6) 
Jo 

where an analytic expression for the integral kernel V A (x,x') is given in 
Eq. Q2] below. This exact result gives the net current as the sum of a baseline 
stationary contribution J A and an excess or "pumped" contribution J ex (x, t) 
associated with the variation of external parameters. This decomposition is 
analogous to that of heat flow encountered in steady-state thermodynamics 



1 From Eq. [2] it follows that J A does not depend on x: in the stationary state, 
the same current flows past every observation point x. 
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[I4lll0|. As we will show, Eq. [6] is amenable to further analysis, and provides 
a useful theoretical tool for studying the response of continuous stochastic 
pumps under arbitrary driving conditions. 

Note that Va(x,x') is not defined uniquely: Eq. [5] is unaffected by a 
transformation of the form Va{x, x') — * Va{x, x') + f(x), since probability is 
conserved. 

To derive Eq. [6l we first solve for P(x, t) in terms of P(x, t) (Eq. [9]), and 
then combine that result with Eq. [3] to determine J(x, t). To this end, let us 
take the following atypical view of the Fokker-Planck equation: for fixed t let 
us interpret Eq. [T]as an operator equation for P(x,t) with operator Ca(x) 
and source term P(x, t). Ordinarily, we solve an operator equation by finding 
the inverse operator. However, since our Fokker-Planck operator has a null 
eigenvector, CaPa = 0> ^ is not invertible. Therefore, we instead introduce 
the integral operator 1Za(x) — J dx'gA(x, x') which is the pseudoinverse of 
Ca{x) (see below for a brief definition of pseudoinverse in this context). Here, 
the integral kernel gA{x,x') is the modified Green's function for Ca{x) [22j . 
defined as the solution of the boundary value problem 



£ A (x')gA(x,x') = 5(x'-x)-P A (x) 



9a(x,x') 



--L 



= 0, 



dg A (x,x') 



-L , (7) 



= 

=0 



dx' 

where C A {x) — Aa{x)8 / dx + BA{x)d 2 / dx 2 is the formal adjoint of C\(x). 
The term P A (x) in Eq. [7] accounts for the fact that £a(x) is not invertible; 
without this term the boundary value problem has no solution [22] . Since 
Eq.[7]is unaffected by a replacement g\(x, x') — > g\[x, x')+f(x), the solution 
of Eq. [7| is not unique. This is the source of the non-uniqueness of V\(x, x') 
mentioned above. 

As mentioned, TZa is the pseudoinverse of Ca- That is, in place of the 
usual inverse property {TZa^-a = 1), TZa satisfies the inverse-like property 

-R A C A P{x,t) = [ dx'gA(x,x')jC. A (x')P(x',t) 

Jo (x) 

= P(x,t) - Pl(x), 

where we have twice integrated by parts and exploited Eq. [7] We see that 
TZaC-a projects onto the orthogonal complement of the null space of Ca 0]- 
Simply put, TZa acts as an inverse on the subspace where Ca is invertible. 

We now apply TZa to both sides of Eq. [T] then use the pseudoinverse 
property (Eq. [5]) to obtain 

P(x, t) = P s A {x) + [ dx'g A {x, x')P(x', t). (9) 
Jo 

Next we apply the current operator (Eq. [3]) to both sides of this equation. 
This gives us 

J(x,t) = J s A + [ dx'JA(x)gA(x,x')P(x',t). (10) 
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Comparing with Eq. [6] we see that Va = Ja3a- Finally, we apply Ja(x) to 
Eq. [7] to arrive at 



C a (x')Va(x,x') = J A (x)S(x' -x)-J s a 



x —L 



o, 



dV A (x,x') 



-l ■ (11) 



x'=0 



x'=o dx' 
This boundary value problems is solved by (see Appendix KJ) 

Va{x,x') = (1 + J A r A (L)) ka{x') + 9(x' -x) + J s a ta(x'), (12) 

where 0{x' — x) is the Heaviside step function; tfA and tpA are given in Eq. 
[SJ and we have introduced the splitting probability 

*a(«) = (13) 
and the conditional mean first exit time 

t a (x)= f dy f V dze My) - VA(z) . (14) 



o Jo 

The splitting probability and the conditional mean first exit time have the 
following interpretations [5]: with A fixed, if the system evolves from xq until 
it first exits the domain [0,-L], then tta{xq) is the probability that this exit 
will occur at x = 0, rather than x = L; and ta(xq) is the average time for the 
system to make this first exit through x = 0. Roughly speaking, the splitting 
probability measures the relative likelihood for the process to go clockwise 
versus counterclockwise around the circle. 

Equations [61 and [T2l specify the current J(x,t), in terms of the rate of 
change of the probability distribution P(x',t). We now investigate conse- 
quences of this result. 

Let us define the excess integrated current <P ex (x) to be the net current 
pumped across the point x, in excess of the time-integrated, baseline station- 
ary flow, <P S = J dtJ s A{ty By Eq.[6l 

& ex (x)= [ dt [ dx'VA(x,x')P(x',t). (15) 
Jo Jo 

If A is varied very slowly from a to b, the system remains near the stationary 
distribution, P(x,t) ~ Pam(x). This suggests that in the adiabatic limit we 

can make the substitution P(x, t)dt — * Va-P a (;e) ■ dA0 in Eq. Qj)]to find 

^(x) = [ dA-E A (a:), (16) 

J a 

where S A (x) = / dx'V A (a;, x')V A P A (x'). This result is geometric: time no 
longer explicitly appears, and the excess current is determined solely by the 

2 This substitution can be justified by appealing to an adiabatic perturbation 
theory. J. Horowitz, S. Vaikuntanathan, S. Rahav, and C. Jarzynski (unpublished). 
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path taken from a to b in parameter space. If the drift and diffusion coeffi- 
cients take the special form A\{x) = (d/dx)VA(x) and Ba(x) = D, Eg. [TBI 
reduces to a result obtained by Parrondo [15]. Analogous geometric expres- 
sions for the adiabatic integrated current have also been derived for discrete 
stochastic pumps [T6 l fT9 j l20 l [2] . 

In Ref. [16] a "no-pumping" theorem was derived for discrete stochastic 
pumps. Briefly, this result stated that in order to obtain non-zero integrated 
currents for cyclic processes, it is necessary to vary both the energy lev- 
els of the discrete states and the effective energy barriers associated with 
transitions between these states. Subsequently, Chernyak and Sinitsyn have 
pointed to the possibility that a similar result applies to continuous stochastic 
pumps [BJ . Here we explicitly develop a no-pumping theorem for continuous 
stochastic pumps. 

Within the general model analyzed above, let us now restrict ourselves 
to the case that detailed balance holds for all A, hence J A = 0. Imagine 
that the parameters are varied periodically with time, from the distant past, 
with period r, so that by t = the system has settled into a periodic steady 
state: P(x, t) = P(x, t + r) 0. In this setup, which is regularly encountered 
in the theory of Brownian motors [3,1,9,17 , the current J(x,t) at a fixed 
location x evolves periodically with time. It is then natural to consider the 
integrated current over one cycle, which for a cyclic process with detailed 
balance ( J A = 0) is 



see Eqs.[3J[ni andlT21 The value of # represents the net circulation of prob- 
ability during one cycle. If the probability merely sloshes back and forth 
without any accumulation of current, then <P = 0; whereas, if # > < 0) 
then there is a nonzero flow of probability in the counterclockwise (clockwise) 
direction. 

We now argue that to obtain <P ^ both the potential ipj\.(x) and the 
splitting probability tt\(x) must be varied during the process. The first of these 
conditions is easy to understand: if the potential remains fixed during the 
process, i.e. (fA(t)( x ) = }5a(s), then the system simply remains in the initial 
equilibrium state, P(x,t) oc exp[— y> a (x)], producing no currents whatsoever; 
this is the "no-go theorem" of Ref. [17] . section 6.4.1. To see that the splitting 
probability must also be varied to produce integrated current, suppose we fix 
7TA(t)(aO = 7r a (a;) but vary ip A ^(x). Then 



since the process is cyclic. We can construct a heuristic interpretation of 
this result by recalling that tta(x) measures the likelihood to generate clock- 
wise rather than counterclockwise flow, as discussed following Eq. 1131 The 

3 As in Ref. [T5 s , our analysis also applies if we begin in the equilibrium distri- 
bution P A ^, then vary A at arbitrary rate around a closed loop in parameter space, 
then allow the system to relax back to equilibrium. 




(17) 




(18) 



7 



integrand n\(x')P(x' ,t) appearing in Eq. [IT] then represents, roughly, the 
contribution to clockwise current induced by the redistribution of probabil- 
ity that occurs at location x' and time t, and the integrated current <P is a 
sum of such contributions. For a cyclic process, any probability that leaves 
the location x' must eventually return, thus if ita(x') remains constant the 
clockwise and counterclockwise contributions ultimately cancel one another 
(<P = 0). If the splitting probability varies with time, then there is no reason 
to expect such cancellation. 

This no-pumping theorem provides a concrete mathematical criterion for 
the generation of zero integrated current. Actually realizing the independent 
variation of both the potential and the splitting probability in any partic- 
ular system will greatly depend on the system's properties and may not 
be feasible; yet, in systems where one may locally vary the drift and diffu- 
sion coefficients independently, one may also vary the potential and splitting 
probability independently, as can be seen from Eqs. [51 and [T3l 

From Eqs.[5land ll31 we see that 7Ta(*) ( x ) [° r equivalently ipA(t) i x )] remains 
constant when the ratio of the drift to the diffusion coefficient does not 
depend on A: 

7TTT = s &- W 

B A (x) 

Thus our no-pumping theorem states that <P = if cither (i) the potential is 
held fixed, or (ii) the drift and diffusion coefficients are related by Eq. [TU 

Finally, we show that Eq. [6] reproduces the known fact that the rectifi- 
cation of current in a periodically driven Brownian motor requires broken 
spatial symmetry [171 15]. Specifically, we will show that when the driving 
protocol is time-periodic and the drift and diffusion coefficients have specific 
spatial symmetries, the integrated current over one period of driving is zero. 
We will say that a periodic function f(x) is symmetric or has even symme- 
try if f(S + x) = f(S — x), or has odd symmetry if f(S + x) = —f{6 — x), 
for some fixed value S. Without loss of generality we set S = 0, since by a 
suitable coordinate shift 6 can take any value. We now assume that the drift 
coefficient has odd symmetry, A\(x) = — A\(— x), and the diffusion coeffi- 
cient has even symmetry, B\{x) — B^—x), for every A. These assumptions 
imply that the system satisfies detailed balance and that both tp\ (x) and the 
periodic steady state P(x,t) = P(x,t + r) are symmetric. Equations [T"3l and 
[TTlthen give 

r T r L r L p^A(y) 

*=/ dt dx dy6{y-x)— L —P(x,t). (20) 

Jo Jo Jo J Q dye^te) 

Changing variables y — > L — y and x — > L — x, and exploiting the symmetry 
and periodicity of ip\(y) and P(x,t), we find 

<P= / dt dx dyeix-y)^ — P(x,t). (21) 

Jo Jo Jo J Q dye^^ 

If we now use the identity Q(x — y) = 1 — 9(y — x) and invoke conservation of 
normalization, j Q dxP(x,t) = 0, Eq. [21] becomes <P = — The anticipated 
conclusion <P = is then obvious. 
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To summarize, Eq. [6] provides a decomposition of the current into a sta- 
tionary contribution, and an excess, "pumped" contribution associated with 
the variation of external parameters. Using this exact result we have estab- 
lished necessary conditions for current generation in cyclic processes - the 
potential tp\(x) and the splitting probability tt\(x) must be varied; demon- 
strated that adiabatic pumping is a geometric effect (Eq. [16]); and verified 
that current rectification requires broken spatial symmetry. Ref. [6] suggests 
the possibility of extending these results to higher dimensions, but this re- 
mains to be done. It would also be interesting to investigate whether the 
no-pumping theorem can be extended to systems without detailed balance. 
We expect that the current decomposition formula will provide a useful tool 
for gaining a deeper understanding of the behavior of stochastic pumps. 
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A Appendix 



The solution to Eq.[Tl]is obtained by combining the homogeneous solution with the 
inhomogeneous solution and then applying the boundary conditions. The two ho- 
mogeneous sol utio ns are the splitting probability which is defined by the boundary 
value problem |24] 







£ A (x')n A (x') 

7T A (0) = 1, 7T A (L) = 



(22) 



and any arbitrary function of x alone, say f(x). The two contributions to the 
inhomogeneous solution, J%ta{x') and 9(x' — x), are obtained by noting that the 
defining boundary value problem for the conditional mean first exit time is |24j 



C A ( X ') ta(x') = -1 
n- (n\ n dr A {L) 

MO) = 0, -gj- 



and that 



t\(x')e(x' - x) 



Aa(x')£j + B A (x') 



dx 

Aa(x') + Ba(x') 



dx'' 



6(x -x) 



JL 
dx' 



Aa{x) 



d_ 

Ox 



B A (3 



= Ja{x)S(x' - x). 
Thus, the most general solution to Eq. II II is 

Va(x, x) = Ctta(x') + f(x) + 0{. 



8(x — x) 
8(x — x) 

) + JXta{x'), 



(23) 

(24) 

(25) 

(26) 
(27) 

(28) 



where C is an arbitrary cons tant. The value of C is fixed by satisfying the first 
boundary condition in Eg. 1111 The second boundary condition is then automatically 
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satisfied due to the structure of Eq. 1111 Lastly, we arrive at the solution in Eq. [12] 
by setting f(x) = 0, which we are free to do since the solution to Eq. [11] is not 
unique. 

To complete the solution we must solve Eqs. [22] and 1231 Both equations can 
be made integrable b y m ultiplying them by the integrating factor e~ VA( - x \ For 
example, consider Eq. [22] for the splitting probability: 

e -v*{*) tf x {x)K A {x) = (29) 

e _w(,) [Aa(x)&; + BaW|t] tta(x-) = (30) 

[_»^«i e -*A(«)^ + e -*A<">£] nA{x) = o (31) 

!^ W ^a(x) = 0, (32) 

whe re t o get Eq.[3l]we used the definition s of y?A and if) a in Eq.[5] By integrating 
Eq. [35] twice we fi nd the solution in Eq. 1131 A similar calculation leads to the 
solution in Eq. [14] for the conditional mean first exit time. 
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